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Abstract. Let A be a random spherical triangle (meaning that vertices 
are independent and uniform on the unit sphere). A closed- form expression 
for the area density of A has been known since 1867; a complicated integral 
expression for the perimeter density was found in 1994. Does there exist a 
closed-form expression for the latter? We attempt to answer this question from 
several directions. An outcome of our work is the exact value of the perimeter 
density at the point ir. 

A spherical triangle A is a region enclosed by three great circles on the unit sphere; 
a great circle is a circle whose center is at the origin. The sides of A are arcs of great 
circles and have length a, b, c. Each of these is < tt. The angle a opposite side a is 
the dihedral angle between the two planes passing through the origin and determined 
by arcs b, c. The angles ft, 7 opposite sides b, c are similarly defined. Each of these 
is < 71 too. 

Define a primal triangle to be a random spherical triangle, obtained by selecting 
three independent uniformly distributed points A, B, C on the sphere to be vertices. 
Define a dual triangle to be a random spherical triangle, obtained by selecting 
three independent uniformly distributed great circles on the sphere to be sides. More 
precisely, starting with independent uniform points A', B', C on the sphere, a dual 
triangle has vertices 

B' xC n A' x C „ A' x B' 
~\\B'xC'\\' ~p'xC"||' ~\\A'xB'\\' 

Hence, while its vertices are not independent, the poles of a dual triangle are. 

Let A be a primal triangle. The area a = a + ft + j — 7r of A satisfies < a < 2tt. 
The perimeter t = a + b + c of A satisfies < r < 2tt. Expressions for the trivariate 
density of (a, (3, 7) and the trivariate density of (a, b, c) are known [TJ but do not give 
useful insight into the distributions of a and r. Crofton & Exhumatus [2] determined 
the density for a: 

(x 2 - Attx + 37r 2 - 6) cos(x) - 6(x - 2ir) sin(x) - 2(x 2 - Airx + 3tc 2 + 3) 

16vrcos(a;/2) 4 
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for < x < 2tt; this formula remained obscure until it was cited in a recent paper 
[3]. Unpublished work of J. N. Boots (mentioned in [1]) was also apparently relevant. 
Jones & Benyon-Tinker [H [5] determined the density for r: 



are complete elliptic integrals of the first and second kind. We wonder: does there 
exist a closed-form expression for this latter density? A direct evaluation of the 
integral does not seem possible, yet conceivably a different geometric argument might 
yield a more accessible formula. We attempt to answer this question from several 
directions. Also, it is clear that the density is zero at x = and diverges to infinity 
at x = 2n. Numerically the density ~ 3\/2/32 to high precision at x = ir, but a proof 
via the preceding is not known. An outcome of our work is a new formula that gives 
the exact value as predicted. 

1. Two Coordinate Systems 
We define two coordinate systems on the unit sphere that will help in our study of 
triangular area and density. 

1.1. Primal Coordinates. Without loss of generality, let A = (1,0,0) and B = 
(cos(k), sin(/t), 0) in xyz coordinates. We wish to locate the unique point C in the 
upper hemisphere so that the triangle ABC satisfies a = 9, b = p, c = k. See Figure 
1. The parameters p, 9 are regarded as varying while the parameter k is fixed. Think 
of rotating the equatorial disk in space so that the vector (1,0,0) remains fixed and 
the vector (0, 1, 0) moves toward (0, 0, 1) through the angle 9. The rotation matrix 
performing this motion is j6] 




for < x < 27T, where 




R 




1 
cos(#) -sin(0) 
sin(0) cos(fl) 
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and 

/ cos(p) \ / cos(p) 
R sin(p) = sin(p) cos(9) 

y J \ sin(p) sin(#) 

which gives the point C. The three-dimensional transformation 



r cos(p) 
r sin(p) 008(6*) 
r sin(p) sin(^) 













- 















has Jacobian determinant 



d(x,y,z) 



d(r, p,6) 



cos(p) — rsin(p) 

sin(p) cos(6 l ) r cos(p) cos(^) — r sin(p) sin(6 ) ) 
sin(p) sin (6*) r cos(p) sin (6*) r sin(p) cos(6 l ) 



r 2 sin(p) 



which implies that the area element in primal p9 coordinates is sin(p)dp d6. 

1.2. Dual Coordinates. Without loss of generality, let A = (1,0,0) and B = 
(cos(p), sin(p), 0) in xyz coordinates. It seems (at first glance) that we should locate 
the unique point C in the upper hemisphere so that the triangle ABC satisfies a — k, 
(3 = 6, c = p. See Figure 2. The parameters p, 9 are regarded as varying while the 
parameter k is fixed. 

Let us examine the great circle containing A, C . It must also contain the point 
(0, cos(k), sin(/c)), since this is the image of (0,1,0) after rotation through angle k. 
Hence a normal vector is V = (1, 0, 0) x (0, cos(k), sin(/«)) = (0, — sin(re), cos(k)). 

Let us examine the great circle containing B, C . Think of rotating the equato- 
rial disk in space so that the vector (cos(p), sin(p), 0) remains fixed and the vector 
(sin(p), — cos(p),0) moves toward (0,0, 1) through the angle 9. The rotation matrix 
performing this motion is j6] 



S 



and 



( cos(p) 2 + (1 — cos(p) 2 ) cos{9) cos(p) sin(p)(l — cos(6')) — sin(p) sin(6 ) ) 
cos(p) sin(p)(l — cos(6')) sin(p) 2 + (1 — sin(p) 2 ) cos(6 ) ) cos(p) sin(^) 
sin(p) sin(^) — cos(p) sin(6 l ) cos(^) 



For example, if p 



sin(p) cos(^) 
- cos(p) cos(6 l ) 
sin(0) 

7r/2, the image of (1,0,0) after rotation through angle 9 is 




(cos(9), 0, sin(^)). As another example, if p = 0, the image of (0, —1, 0) after rotation 
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through angle 9 is (0, — cos(6 l ), sin(#)). Hence the great circle must contain the point 
(sin(p) cos(6>), — cos(p) cos(6 l ), sin(6 1 )) and a normal vector is 



W 




sin(p) cos(#) 
x | — cos(p) cos(6) 
sin(0) 



sin(p) sm(8) 
— cos(p) sin(6 l ) 
- cos(0) 













- 















The point C is orthogonal to the two normal vectors and at unit distance from the 
origin, equivalently, C = (V x W)/ \\V x W\\. 

Now, in fact, this is more than what is required. We need only (on second glance) 
specify the great circle containing B, C and this is done via locating W or — W. The 
three-dimensional transformation 

— r sin(p) sin(#) 
r cos(p) sin^) 
r cos(#) 

has Jacobian determinant r 2 sin(#), which implies that the area element in dual p9 
coordinates is sin(8)dpd9. Perhaps this is obvious by duality. It is good, however, to 
see the supporting geometric details. 

2. Four Approaches 
We illustrate using four different trigonometric identities and the above two coordi- 
nate systems. More possible approaches will be mentioned in a later section. 

2.1. Primal Area. As in section [l.lj . assume that the triangle ABC satisfies 
a = 9, b = p, c = k. These three parameters are related to area a as follows: 



tan 



cot 



esc ( 9 



sin 



A proof appears in section [6.1] . For fixed a and k, define 



/(*) 



7T 

2 arctan 



cot 



esc ( 9 



a 



sin 



if < 9 < a/2, 
if a/2 < 9 < n 



then the conditional probability, given c, is 

j TT f(8) 

P {area < a \c = k} = — 

2n 



sin(p)dp d9 

— J J sin (p) dp d9 + — J J sin(p)dpd9 

a/2 





a/2 7^ 
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1 - cos(/(0)))d0j . 

This result can be experimentally verified by generating many primal triangles ABC 
with c = k, and then plotting all pairs (6, p) corresponding to triangles with area 
< a. The scatterplot fills the region [0, 7r] x [0, 7r] except for the portion lying above 
the curve p = f(9). 

We will later discuss |3.1] how the unconditional probability P {area < a} is eval- 
uated exactly, for arbitrary a. The method is quite long and intricate. 

Here is a quick method for evaluating not probability, but instead density, at 
a = 7i. We start with the conditional density 



d 1 
da 2n 



a 



1 - cos(f(6)))d8 



)-*(■ 



a/2 



1-cos / § 



— (1 - <xx(f(0)))M 
da 



a/2 



— J sm{f{9))g{6)d9 

a/2 



where 



fdn(/(0)) 



2 tan (f) sin (O - f) sin 
tan (f) 2 sin (0 - f) 2 + sin(f 



since 



sin(2 arctan(C)) 



2C 



and where 



4f 



l + C 2 

tan (j) sin(0) 



rf(J tan(f) sin(0-f) +sin(f)' 
It follows that the unconditional density is 



I j J tan (f ) sin (g - |) sin (f ) sin (6) . 

' ' — : — 2 sin(/cjd0 are 



2tt 



CT /2 



tan(f) sin(0-|) + sin (f 



because the density for k is sin(fc)/2. By the half-angle formula for tangent, this is 
the same as 



1 

2~it 



a/2 



1 - cos(k)) (1 + cos(k)) sin (6 - f ) sin (f ) sin (6) 
— — 2 sin(re)(i0 dn. 



(1 - cos(«)) sin (0 - |) 2 + (1 + cos(«)) sin (|)" 
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In the special case when a = tt, this becomes 



1 

Ytt 



cos (9) sin (9) (1 — cos(k) 2 ) sin(/«) 



tt/2 

1 



1 + cos(k) + (1 — cos(k)) cos (Qy 



-A8 dn 



2tt 



u (1 — V 1 



i x i x [l + v + (l-v) u 2 ] 



irdu dv = — 
2 Att 



consistent with Crofton & Exhumatus. No analogous simplication seems to occur, 
for example, when a = n/2 or a = 37r/2. 

2.2. Dual Perimeter. As in section [1.2] , assume that the triangle ABC satisfies 
a = k, j3 = 8, c = p. These three parameters are related to perimeter r as follows: 

tan r-) 



cot I — sin p esc — I . 



A proof appears in section [6.3] . For fixed r and k, define 

cot (f) sin(|-p) esc (|) 



f(p) 



2 arctan 




if < p < r/2, 
if r/2 < p < tt 



then the conditional probability, given a, is 

TT f( P ) 



P {perimeter < r \a = n} 



1 

2^ 



sm(8)d8dp 

-J J sin(6)d6 dp + — J Jsm(8)dddp 
oo 

r/2 

— f (1 - cos(/(p)))dp. 





r/2 /(p) 



7T 



r/2 



This result can be experimentally verified by generating many dual triangles ABC 
with a = k, and then plotting all pairs (p, 8) corresponding to triangles with perimeter 
< r. The scatterplot fills the region [0, tt] x [0, tt] except for the portion lying above 
the curve 9 = f(p). 

Since (dual perimeter) = (27r— primal area), it is not surprising that conditional 
probabilities are so similar. 
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For completeness' sake, let us compute the conditional density 



d 1 
dr 2tt 



t/2 



1 - cos(f(p)))dp 



1 
1 

2^ 



1 - cos (/ (J) 



t/2 



+ 



\ 



1 - cos(f(p)))dp 



I 



t/2 



where 



2 tan f f J sin ( | - p) sin 
tan f |) sin f |) + sin (| — p) 

tan (|) sin(p) 
tan (|) sin + sin (| — p) 
It follows that the unconditional density is 



sin(/(p)) 



f \ d f 
9{P) = Tt 



2 ■ 



1 

2^ 
1 

27 



TT T/2 



tan f|) sin (| — p) sin sin (p) 



o o 

TT T/2 



sin(n)dpdK 



tan(f) sin(j) + sin (| - p) 
1 — cos(/t)) (1 + cos(k)) sin (| — pj sin sin (p) 



o o 



sin{n)dpdK 



(1 — cos(k)) sin + (1 + cos(k)) sin (| — p) 
because the density for k is sin(fi;)/2. In the special case when r = tt, this becomes 



1 

2^r 



tt vr/2 



cos (p) sin (p) (1 — cos(k) 2 ) sin(fi;) 



1 

27 



o 

1 i 



1 — cos(k) + (1 + cos(k)) cos (pY 



-dpdK 



u(l-v 2 ) 



-1 



\l-v + (l + v) V?? 



du dv 



An 



consistent with Crofton & Exhumatus. 



2.3. Primal Perimeter. As in section jl.lj . assume that the triangle ABC sat- 
isfies a = 8, b = p, c = k. These three parameters are related to perimeter r as 
follows: 

sin(r — k) cos(r — k) — cos(k) . , 
cos(#) = — , ; + i — ^ — cot(p). 



sin k 



sin k 
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A proof appears in section |6.2] . For fixed r and k, define 
f(p) = 




sin(r— k) . cos(r— k)— cos(k) 
sin(re) sin(fc) 



if 0<p<r/2-/t, 
cot(p)] if r/2 - k < p < r/2, 
if r/2 < p < 7r 



assuming k < r/2; otherwise f(p) = 0. Then the conditional probability is 

TT f(p) 



P {perimeter < r |c = k} 



1 
2^ 

1 

2tt 
1 

2^ 



sin(p)<i# dp 



o o 

( T/2-K TT 



r/2 f(p) 



TT 



\ 



sin (p) cf# dp + y y sm(p)d9dp+ J J sin(p)d9dp 



r/2-K o 
r/2 



r/2 



7T 



V 



1 — cos ( £ — 



(5-*) 



+ 



/(p) sin(p)dp 



t/2-k 



This result can be experimentally verified by generating many primal triangles ABC 
with c = k, and then plotting all pairs (p, 9) corresponding to triangles with perimeter 
< t. The scatterplot fills the region [0, 7r] x [0, n] except for the portion lying above 
the curve 9 = f(p). 

An exact evaluation of the unconditional probability P {perimeter < r}, for arbi- 
trary r, remains open |3.2j . 

Here is a quick method for evaluating not probability, but instead density, at 
t = 7i. We start with the conditional density 



d 1 
dr2n 



( 



r/2 



f(p) sin (p) dp 



7-/2-K 



1 — COS f | — k) 
± (f sin (!-«) + |/ (l) sin (|) - i/ (| - «) sin (f - K 



r/2 



-^-/(p) sin(p)dp 



r/2-/* 



1 

2^ 



r/2 



#(p) sin (p)dp 



T/2-K 



where 



M = -jz 



sin(r — k — p) sin(p) 



^ r \/ sin(fi;) 2 sin(p) 2 — [cos(«) cos(p) — cos(r — k — p)\ 
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It follows that the unconditional density is 

1 sin(r-ft-p)sin(p) ..^ K ) dpdK 



^ o t/2-k V srn ( K ) 2 sin(p) 2 — [cos(ft) cos(p) — cos(r — k — p)\ 
because the density for k is sin(/c)/2. In the special case when r = n, this becomes 

tt/2 tt/2 

Iff sin(ft + p) sin (ft) sin(p) 

-.dp an 



lrc { n/ J 2 _ K yj sin(ft) 2 sin(p) 2 - [cos(ft) cos(p) + cos(ft + p)f 
Iff sin (ft + p) smiK) sm(p) 

dp dK 



4tt 



W2- K V~ 4 COS ^^ COS( ^^ C ° S ^ + P ) 



vr/2 tt/2 

Iff sin(ft + p) sin(ft) sin(p) 

-dpdK 



>7T 



tt/2-k 



— cos(ft + p) J cos(ft) W cos(p) 



Iff sin(ft + p) sin(ft) sin(p) 

dpdn 



47T 



7r/4 7r /2—k \f~ C °^ K + S ^/ C0S ( K ) COS (P) 

Iff sin(ft + p) cos(ft — p) — cos(ft + p) 

ap aft. 



^^J/^/Lk ^/- cos(ft + p) v /cos(k-p) + cos(ft + p) 

Let u — K + p, v — k — p. Then v)/8(k, p)| = 2 and the integral is transformed 
to 

7r/2 tt—V . . , — 

sin(u) cos(t>) — cos(w) , , 3\/2 

era ctf 



8^7? V J ^— cos(-u) ycos(f) + cos(w) 32 



as promised. 



2.4. Dual Area. As in section |1.2j . assume that the triangle ABC satisfies a = 
ft, (3 = 6, c = p. These three parameters are related to area a as follows: 

. , sinfcr — ft) cosfcr — ft) — cos(ft) ... 

- cos p = \ > + —± rf- ^ cot . 

sin (ft) sin (ft) 

A proof appears in section [6.4j . For fixed a and ft, define 

7T if0<#<<x/2, 
f(ti) = { vr - arccos + cot(0)] if a/2 < 9 < n - (ft - a/2), 

if 7T- (ft -a/2) < 9 < IT 
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assuming k > a/2; otherwise f{6) = ir. Then the conditional probability is 



7T f(0) 



P {area < a \a = k} = — J J sin(6)dp 



dO 



1 

2~ir 
1 

2^ 



o o 

(a/2 



7r-(«-<r/2) f(8) 



TV 



sin(9)dpd9 + 



sm(6)dpd6 + 



sin(9)dp d6 



\ o o 



7T 



1 - 008 (f)] + 



a/2 
7r-(«-<r/2) 



n-(K-a/2) 



\ 



/(0) sin(0)d0 



cr/2 



This result can be experimentally verified by generating many dual triangles ABC 
with a = k, and then plotting all pairs (0, p) corresponding to triangles with area 
< a. The scatterplot fills the region [0, it] x [0, it] except for the portion lying above 
the curve p = f(8). 

Since (dual area) = (2n— primal perimeter), it is not surprising that conditional 
probabilities are so similar. 

For completeness' sake, let us compute the conditional density 



d 1 
da 2tt 



■K-(n-a/2) \ 

f{9) sin(0)d0 

J 



a/2 



1 - COS (f ) 

= i (f sin (f ) +1/(^-^-1)) sin - (« - f )) - 5/ (f ) sin (f 



w-(k-<t/2) 



a/2 

7T-(«-<7/2) 



(i<7 



/(0) sin(0)d0 



^ | <?(0)sin(0)d0 



where 



cr/2 



sin(o" — k — 6) sin(0) 



^ ^sin(«;) 2 sin(0) 2 — [cos(k) cos(0) — cos(cr — k — 9)\ 
It follows that the unconditional density is 



1 

Air 



7T 7r-(ft-(T/2) 



sin((T — k — 9) sin(0) 



cr/2 cr/2 



sin(/t) 2 sin(0) 2 — [cos(/c) cos(0) — cos(o" — /c 



sin(/t)(i6 l dft 
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because the density for k is sin(ft)/2. In the special case when a = it, this becomes 

7T 37T /2 — K 

lrr sin(ft + 9) sin(ft) sin(6>) 
/ / ^=d0 an 

J/2 J/2 \/sin(/t) 2 sin(#) 2 - [cos(ft) cos(#) + cos(ft + 9)f 

7T 37r/2 — K 

Iff sin(ft + #) sin(ft) sin(6 l ) 

-du du 



47V Jo Jo a/ -4 COS ( K ) COs(#)cOs(ft + 



tt/2 tt/2 
7r 3-7r/2— re 

Iff sin(/t + 0) sin (/c) sin(0) 
— / / — ^^^^=^^^^= — = — =do aft 

^2 tt/2 V /- COS ( K + ^) V-COS(K) ^-COS(0) 

3tt/4 3tt/2-k 

Iff sin(ft + 0) sin(ft) sin(#) 

-du aft 



47T 

tt/2 



■[ ^ — cos(ft + 9) y / cos(ft) cos(#) 



3tt/4 3tt/2-k 

I /■ /■ sinfft + 9) cosfft — 0) — cosfft + 9) 

=d9 aft 



4 v / 2 ?r J /2 { y / -cos(ft + 6») y / cos(ft - 0) + cos(ft + 6 

7r/2 37r/2 _ 

1 /" /" sin(w) cos(v) — cos(m) , , 3v2 

' dudv - 



J q yC C0S ( M ) y/ cos(f ) + COs(u) 32 

as promised. 

3. Two Evaluations 
3.1. Successful Evaluation for Primal Area. Starting from the half-angle 
formula for tangent 

ta ,pV l-cos(p) 



2/ l + cos(p) 
we deduce that 

cot m' + 1= l±£^4 + 1 _ 2 



,2/ 1 — cos(p) 1 — cos(p) 

hence 

1 — cos(p) 1 1 



2 cot(f) +1 fi 2 sin(#-f) +1 

by |2.1j . where 



= tan [ — ) esc [ — 1 . 
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It follows that 

" 1 - COs(/(0)) 



dd 



a/2 



a/2 



n 2 sin(^-fj +1 



-dO 



— arctan 

arctan ( \/Vt 2 + 1 tan(cr/2) 



V^ 2 + l 



if < a < 7r, 

if 7T < (T < 27T 



and thus 



P {area < a} 



1 

2^ 
1 

2^ 



7T 






o 


/ 


2 







7T 






(T 


/ 


2 








7T 



— arctan ( + ltan(cx/2) 



sm(K,)dK if < a < 7r, 
sin(/c)dK if 7r < a < 27r. 



The area density is therefore 



1 

2^ 
1 

2^ 



^ 7 7r — arctan ( \fVt 2 + 1 tan(cr/2 

H / v , 

da J JQ 2 + 1 



d J arctan (V^ 2 + 1 tan(<x /2) 



(icr 



sin(K,)dK 
sin(/«)c?K 



if < (X < 7T, 



if 7T < a < 27T 



which, as outlined in |7.1j . gives rise to the Crofton/Exhumatus expression. 

3.2. Unsuccessful Evaluation for Primal Perimeter. There does not seem 
to be an analogous approach for computing 



r/2 



f(p) sin (p)rfp 

t/Z-k 

from |2.3] in closed-form. We suspect that elliptic integrals will arise, but have not 
yet found a method for demonstrating this. See |7.2] for more details. 

4. Two More Approaches 
4.1. Median Area. Assume that the triangle ABC satisfies c = k. A median 

in ABC is the great circle drawn from vertex C to the midpoint P of side c. Let p 
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denote the spherical distance between P and C, and 6 denote the angle between PB 
and PC. These three parameters are related to primal area a as follows: 

sin(/c/2) sin(p) sin(6 l ) 



a 

tan — 



cos(k/2) + cos(p) 
A proof appears in section [6.5] . For fixed a and k, define 

(f) 



/(«) 



arccos 



cos 



tan 



(f) csc(0) s 



=F sin ~ 



tan 



(f) csc(0)= 



COS 



(f) -tan(f) csc(#) 2 -l 



where — is chosen if a < n and + is chosen if a > n; then the conditional probability 
is 



P {area < o \c = K } 



1 

2^r 
1 

2^ 



f(9) 



sin(p)dpd9 



o o 



1 - cos(f(6)))de 



I j( cos(f)tan(f) 2 csc(^) 2 T-sin(f) 2 yi + tan(f)"csc(^)^ 
27 1 



V 



cos 



(f) -tan(f) csc(#) 2 -l 



dO. 



This result can be experimentally verified by generating many primal triangles ABC 
with c = k, and then plotting all pairs (8, p) corresponding to triangles with area < a. 
The scatterplot fills the region [0, ir] x [0, tt] except for the portion lying above the 
curve p = f(9). This approach is believed to be the same as Crofton & Exhumatus 
(details in [2] are rather thin). Not seeing any advantage over our approach in [2.1] , 
we stop here. 

4.2. Bisector Perimeter. Assume that the triangle ABC satisfies a = k. An 
angle bisector in ABC is the great circle drawn from vertex A that splits angle a 
in half. Define Q to be the intersection point between this circle and side BC. Let p 
denote the spherical distance between Q and A, and 9 denote the angle between QC 
and QA. These three parameters are related to dual perimeter r as follows: 

/t\ cos(k/2) sin(p) sin(^) 

ai1 V2 ) ~ ~sin(fi;/2) + cos(p) sin(0)' 

A proof appears in section |6.6] . For fixed r and k, define 

tan(r/2) sin(/t/2) 



flower (p) 



arcsm 



tan(r/2) cos(p) + cos(k/2) sin(p) 
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/upper (p) = 7T - arcsin 

assuming 



tan(r/2)sin(/c/2) 



tan(r/2) cos(p) + cos(k/2) sin(p) 



/ cos(t/2) + Bin(«/2) \ ^ ^ 

arccos -— , = p thrcs < p < tt. 

\ 1 + cos(r/2) sm(K/2j / 

The region of all pairs (p, 9) corresponding to triangles with perimeter < r is more 
complicated than earlier examples. The scatterplot fills the region [pthres, 71 "] x [0, tt] 
except for portions lying either above the curve 9 = / upper (p) or below the curve 
9 = /lower (p)- The conditional probability is 



P {perimeter < r \a = k,} = — J J sm(9)d9dp 

/lower (p) 

(cos(/i owcr (p)) - cos(/ upper (p)))dp 



1 

2^ 



res 



tau(r/2) sin(«/2) 



dp. 



tt J \ l v tan(r/2) cos(p) + cos(k/2) sin(p) 

Pthrcs 

Due to the unanticipated complexity, we stop here. 

5. Two More Coordinate Systems 
The primal coordinate system jl.lj allows us to specify a triangle, given a fixed side k, 
with an additional side p and an angle 9. Can we do as well with two angles instead? 
The dual coordinate system |1.2j allows us to likewise specify a triangle, given a fixed 
angle k. Can we do as well with two sides instead? 

5.1. Angle Coordinates. Without loss of generality, let A = (1,0,0) and B = 
(cos(k), sin(/c), 0) in xyz coordinates. We wish to locate the unique point C in the 
hemisphere so that the triangle ABC satisfies a — (p, (3 — ip, c — k. See Figure 3. 
The parameters (p, ip are regarded as varying while the parameter k is fixed. 

Let us examine the great circle containing A, C . It must also contain the point 
(0, cos(<p), sin(<p)), since this is the image of (0, 1,0) after rotation through angle (p. 
Hence a normal vector is V = (1, 0, 0) X (0, cos(<p), sin(<p)) = (0, — sin(<p), cos(<p)). 

Let us examine the great circle containing B, C. Think of rotating the equato- 
rial disk in space so that the vector (cos(k), sin(/«), 0) remains fixed and the vector 
(sin(/t), — cos(k), 0) moves toward (0, 0, 1) through the angle ip. The rotation matrix 



Random Spherical Triangles 



15 



performing this motion is [6] 



cos(k) 2 + (1 — cos(k) 2 ) cos('0) cos(k) sin(fi;)(l — cos(ip)) 



sin(/t) sin(ifj) 



S 



cos(k) sin(/t)(l — cos(-?/>)) 
sin(/t) sm(ip) 



and 



sin(ft) 2 + (1 — sin(/t) 2 ) cos(^) cos(k) sm(ip) 

— cos(k) sin(ip) cos(^) 

sin(fi;) cos(V') 
— cos(k) cos( , 0) 
sin(-0) 

Hence the great circle must contain the point (sui(k) cos(ip), — cos(k) cos (ip), sin (?/>)) 
and a normal vector is 





( 


sin(«;) 




s 




— cos(k) 


)-( 




\ 










( cos(k) ^ 




W = 


sin(«) 


•( 




V o J 





sin(/t) cos(^) 
- cos(k) cos(^) 
sin(^) 



/ sin(/«) sin(^) 

— cos(k) sin(^) 
\ — cos(ip) 



The point C is orthogonal to the two normal vectors and at unit distance from the 
origin, equivalently, C = (V x W)/ \\V x W\\. We have 



V x W 



I sin(yj) cos("0) + cos(k) cos(y?) sin(/0) 

sin(/c) cos(<^) sin( , 0) 
y sin(ft) sin (<^) sin(^>) 



\\V x W\ = \fl — (cos(<p) cos(t/>) — cos(k) sin(y?) sin(?/>)) 2 
and thus the Jacobian determinant of (r, (p, ip) 1— > (x, y, z) = —rC simplifies to 

sin(ft) 2 sin(</?) sin(?/>) (sin(</?) cos('0) + cos(k) cos(</?) sin(-?/>)) 2 + sin(fi;) 2 sin^) 2 



1 — (cos((^) cos('0) — cos(k) sin(y) sin(^))" 



5/2 



5.2. Side Coordinates. Without loss of generality, let A = (1, 0, 0) and B = 

(cos(£), sin(£), 0) in xyz coordinates. It seems (at first glance) that we should locate 
the unique point C in the upper hemisphere so that the triangle ABC satisfies a — k, 
c — £, b — 77. See Figure 4. The parameters £, 77 are regarded as varying while the 
parameter k is fixed. Think of rotating the equatorial disk in space so that the vector 
(1, 0, 0) remains fixed and the vector (0, 1, 0) moves toward (0, 0, 1) through the angle 
K. The rotation matrix performing this motion is [6] 



R 
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cos(?7) 
cos(k) sin(Ty) | , 
sin(re) sin(r/) 



and 

/ cos(?7) 
R sin(?7) 

V o ' 

which gives the point C. 

Now, in fact, this is less than what is required. We must (on second glance) 
specify the great circle containing B, C. This is done via a normal vector U = 
(B x C)/\\B x C\\, where 



BxC 



' cos (0 \ / cos(r]) 

sin(£) x cos(k) sin(r/) 
\ / \ sin(K) sin^) 

/ sin(/t) sin(£) sin(7/) 

— sui(k) cos(£) sin(^) 
\ — sin(£) cos(?7) + cos(k) cos(£) sin(r/) 



\B x C|| = y 1 — (cos(£) cos(r/) + cos(k) sin(£) sin(?7)) . 
Thus the Jacobian determinant of (r, £, if) h-> (x, y, z) = rlf simplifies to 

sin(fi:) 2 sin(£) sin(?7) (sin(£) cos(?7) — cos(k) cos(£) sin(^)) 2 + sin(ft) 2 sin(^) 2 



1 — (cos(£) cos(?7) + cos(k) sin(£) sin(r/)) z 



5/2 



5.3. Possible Applications. Combining an identity in [7] with the Law of Cosines 
for Angles, we obtain 



tan 



T 



sin(y?) sin(?/>) sin(/t) 



.27 cos(<^) + cos(^) — cos(v?) cos(-?/>) + cos(k) sin((/9) sin(^) — 1 
Let us solve for ip as follows: 



[cos(y) + cos(^) — cos(ip) cos(^) — 1] tan ( — ) = sin(y?) sin(-?/>) 



T 



sinf/v) — cos(/.') tan I — 



hence 



[1 — cos(y?)] [1 — cos(/0)] tan ( — ) = sin(<£>) sin(-?/>) 



T 



sin(ft) — cos(k) tan 



T 



hence 



1 — cos(y?) / t\ sm(?p) 
— - — t an I — J = — — 

sin(<^) V27 1 — cos(^) 



sin(/.') - cos(/.') tan. ( — 
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hence 



hence 



tan 



esc 



K sin 



sin (■0) 
1 — cos(ip) 



cos 



tan £ csc U-k sin m -1 



t 



an (|) esc (j - ft) sin (j) + 1 



because y = y/l — x 2 / (1 — x) has inverse x = (y 2 — 1) / (y 2 + 1). For fixed r and k, 
define 

tan(f) csc(f-tt) sin(f) -1 



arccos 




tan(^) csc(^— k) sin(-|) +1 



if < k < r/2, 
if r/2 < k < vr 



then the conditional probability, given c, is 



f(e) 



P {perimeter < r |c = k 



iff sin(ft) 2 sin((^) sin(i/>) [(sin(ip) cos(i/j)+cos(k) cos(</p) sin(-i/))) 2 +sin(K) 2 sin(i/>) 2 ] , , 







1— (cos(ip) cos(i/))— cos(k) sin(y) sin(i/)))' ! 



Similarly, combining an identity in |7J with the Law of Cosines for Sides, we obtain 

sin(£) sin(?7) sin(«;) 



tan 



a 



1 + cos(£) + cos (7/) + cos(£) cos(r/) + cos(k) sin(£) sin(^) 



Let us solve for rj as follows: 



[1 + cos(£) + cos(r]) + cos(£) cos(?y)] tan ( — ) = sin(£) sin(ry) 



a 



sin k 



cos( /v } l an | — 



hence 



hence 



1 + cos(£) fa 
— ■ f{: \ tan 



sin(?7) 



siniK) — cos(k) tan 



cos (7/) 



1 + cos (77) 

1 — cot (Jj\ esc (k — |) sin 

1 + cot (|) 2 esc (k - f ) 2 sin (§)' 



because y = yl — x 2 / (1 + x) has inverse x — (1 — y 2 ) / (1 + y 2 ). For fixed a and k, 
define 

vr ifO<K<tr/2, 



/(0 



arccos 



1— cot(^) csc(k— ^- 1 sin 

i +cot (i) 2 csc ( K _!) 2 sin (!) 2 



if a/2 < k < vr 
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then the conditional probability, given a, is 

i n m r 2 i 

, | ..iff sin(K) 2 sin(£) sin(rj) (sin(£) cos(r))— cos(re) cos(§) sin(»?)) +sin(K) 2 sin(»7) 2 

P {area <<J\a = K\ = — I I : — 1 --572 L ar] a£. 

Ill J J l-(cos(£) cos(»j)+cos(k) sin(£) sin(r?)) 

L 1 

We have not further pursued this direction of inquiry. 

6. Trigonometric Identities 
The following formulas are required in the main text. 

6.1. Primal Case i. To prove 



fb\ fc\ f o\ . (a 

tan — = cot - esc a - - sin 



2) \2J V 2/ V2 

we expand cos(cr/2) and make use of Delambre's analogies [S]: 



cos 



.2, 



cos 



j3 + 7 7r — a 



cos — sm y +sm — cos j^- 



(b + c\ . fa 
cos I— sm f- 



2 (b-c\ fa 

+ cos — - — COS 



sec 



/ b c b c\l — cos a ( b c b c\l + cos a 

cos - cos sm - sm h cos - cos — h sm - sm - 

\ 2 2 2 2) 2 \ 2 2 2 2) 2 

b c . b . c \ a 
cos - cos — h sm - sm - cos a sec — . 
2 2 2 2 2 



Also 



sin f — 

V2, 



sm 



'/3 + 7 



7r — a 



+ (a\ . //3 + 7 \ . fa 

~ C ° S [— ) C ° S { 2 J + Sm J Sm U 

fb + c\ fa\ fa\ (b-c\ fa 

-cos 1 — 1 sm i-j cos (-) + cos 1 — 1 cos f - | sn, 



'6c be 

cos - cos sm - sm - 

2 2 2 2 



6 c . 6 . c' 
cos - cos — h sm - sm - 
2 2 2 2 





'a 


sec | 


k2 


a . 


a 


— sm 




2 


2 



a 
2 



6.c a. a a . 6 . c . a 
2 sm - sm - cos — sm — sec — = sm - sm - sm a sec — . 
22222 22 2 
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Dividing, we obtain 



hence 



and therefore 



a cos 77 cos 75 + sin 75 sin 75 cos a 

cot - = \ 

2 sin 2 sin | sin a 



o . (be \ . a 

cos — sin a = cot - cot — h cos a sin — 
2 \ 2 2 / 2 



a 



<j .a b c . a 

an a = sin a cos cos asm - = cot - cot - sin — 

2/ 2 2 2 2 2 



as was to be shown. 

6.2. Primal Case ii. To prove 



. . sin(r — c) cos(r — c) — cos(c) 

cos a = — K — + ^ — cot (6) 

sm(cj sm(cj 

we expand cos(a): 

cos(a) = cos(r — b — c) 

= cos (—b + (r — c)) 

= cos(6) cos(r — c) + sin(6) sin(r — c) 

and make use of the Law of Cosines for Sides: 

cos(a) = cos(6) cos(c) + sin(6) sin(c) cos(a) 

thus 

sin(6) sin(c) cos(o;) = sin(r — c) sin(6) + [cos(r — c) — cos(c)] cos(6). 
Alternatively, we have 

sin(c) cos(a) = sin(r — c) + (cos(r — c) — cos(c)) cot(6) 

hence 

. , cos(r — c) — cos(c) 

-tan(6) 



sin(r — c) — sin(c) cos(a;) 
but solving for b turns out to be more complicated than our strategy of solving for a. 
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6.3. Dual Case i. To prove 




we expand cos(r/2) and make use of Delambre's analogies [8]: 



008 © = cos (^ + i) 

fa + b\ fc\ (a + b\ (c 

= cos I — cosy -sm I— sm(~ 



(a + B\ (c 
cos I— COS I- 



— COS 



a — (3 s 



sm 



/ a 8 a B\ 1 + cos c 

cos — cos sm — sm — 

V 2 2 2 2/2 

. a . 8 a 8 \ 7 
— sm — sm — h cos — cos — cos c esc — . 
2 2 2 2 /2 





) 2 








a 


cos 






2 



CSC 



8 a 3 s 

h sm — sm — 

2 2 2 



cose 



7 

esc — 
2 



Also 



Sin = Sin (^ + ^j 

fa + b\ . fc\ . fa + b\ (c 
= cos I— smf-)+sm I— cos f - 



(a + 8\ (c\ . (c 

cos cos - sm - 

\ 2 \2J V2 



fa-8\ . fc\ (c 

cos sin - cos - 

\ 2 \2J \2 



a 8 . a . 8\ I a 8 . a . 8 

cos — cos sm — sm — + cos — cos — h sm — sm — 

v 22 22/^22 2 2, 

a 8 c . c 7 a 8 . 7 

2 cos — cos — cos - sm - esc — = cos — cos — sm c esc — . 
22222 22 2 



CSC 



c . c 7 
cos - sm - esc — 
2 2 2 



Dividing, we obtain 



r — sin 77 sin § + cos 77 cos £ cos c 

COt - = ? 2 7f J 

2 cos 7j cos | sin c 



hence 



T / a 8 \ . r 

cos — sm c = — tan — tan — h cos c sm — 
2 \ 2 2 / 2 
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and therefore 



. (t \ . t t . a 8 . T 

sin c = sin — cos c — cos — sin c = tan — tan — sin — 

V2 J 2 2 2 2 2 

as was to be shown. 

6.4. Dual Case ii. To prove 

. , sin(<j — a) cosfcr — a) — cos(a) 

- cos c = — ) . , ; + ^ r^- K —L cot /3 

sm(aj sin (a J 

we expand — cos(7): 

— cos(7) = — cos((j — a — 8 + 7r) 

= cos (— 8 + (cr — a)) 

= cos(/3) cos(cr — a) + sin(/3) sin(cr — a) 

and make use of the Law of Cosines for Angles: 

— cos(7) = cos(a) cos(/3) — sin(a) sin(/3) cos(c) 

thus 

— sin(a) sin(/3) cos(c) = sin(cr — a) sin(/3) + [cos(a — a) — cos(a)] cos(/3). 
Alternatively, we have 

— sin(a) cos(c) = sin(<T — a) + (cos(a — a) — cos(a)) cot (8) 

hence 



cos(a — a) — cos(a) 

-tan(/3) 



sin(cr — a) + sin(a) cos(c) 
but solving for /3 turns out to be more complicated than our strategy of solving for c. 

6.5. Median Case. Let p, 6 be defined within triangle ABC as in |4.1) . Applying 
the Law of Cosines for Sides to both triangles CPB and CPA, we have 

cos(a) = cos(p) cos(c/2) + sin(p) sin(c/2) cos(6 l ), (1) 

cos(6) = cos(p) cos(c/2) — sin(p) sin(c/2) cos(9) (2) 
because cos(7r — 6) = — cos(#); hence 

1 + cos(a) + cos(6) + cos(c) = 1 + 2 cos(p) cos(c/2) + cos(c) 

= 2 cos(p) cos(c/2) + 2cos(c/2) 2 
= 2 cos(c/2) (cos(p) + cos(c/2)) . 
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Applying the Law of Sines to both triangles CPB and ABC, we have 

sin(a) sin(p) sin(6) sin(c) 
sin(#) sin( ( 5) ' sin(/3) sin(7) 

hence 

sin(p) sin(0) sin(/3) sin(c) . 

sin a sin (6 sin 7 = — — —— sin 7 

sin(p) sin (7) 

= sin(p) sin(0) sin(c) 

= 2sin(p) sin(0) sin(c/2) cos(c/2). 

Eriksson [7] proved that 

'cr\ sin(a) sin(6) sin(7) 
tan [ — 



. 2 / 1 + cos(a) + cos(6) + cos(c) 
from which 

/er\ sin(c/2) sin(p) sin(0) 
tan V2 J ~ cos(c/2) + cos(p) 

follows immediately. 

Adding equation (j2J) to ([T]), we obtain [!) 

. . cos(a) + cos(6) cos (ef) cos fe) 

2 cos (j J cos(j) 
subtracting equation §2$) from ([!]), we obtain 

... cos(a)-cos(6) sin (^) sin (^) 

COS(6>) = -^r = ' ^ '-. 

2sin ( |J sin(p) sin(|Jsin(p) 

Thus, given a, b, c, it is easy to compute p and then 9 (in that order). 
Rearranging equation (J3J) to 



tan (a/2) sin(c/2) sin(p) y/1 - cos(c/2) 2 yl - cos(p) 2 
sin(0) cos(c/2) + cos(p) cos(c/2) + cos(p) 

that is, 

z = 

y + x 

we solve for x = cos(p): 

(y + x) 2 z 2 = (l-V 2 ) (l-x 2 ), 

that is 

(l - y 2 + z 2 ) x 2 + (2?/2 2 ) x - (l - y 2 - y 2 z 2 ) = 
and obtain the expression for p = f(9). 
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6.6. Angle Bisector Case. Let p, 9 be denned within triangle ABC as in [4.2j . 
Applying the Law of Cosines for Angles to both triangles AQC and AQB, we have 

— cos(7) = cos(#) cos(a/2) + sin(#) sin(a/2) cos(p), (4) 

— cos(/3) = — cos(#) cos(a/2) + sin(0) sin(a/2) cos(p) (5) 
because cos(7r — 9) = — cos(9) and sin(7r — 9) = sin(6 l ); hence 

cos(a) + cos(/3) + cos(7) — 1 = cos(a) — 2 sin(0) sin(a/2) cos(p) — 1 

= -2 sin(0) sin(a/2) cos(p) - 2 sin(a/2) 2 
= -2 sin(a/2) (sin(0) cos(p) + sin(a/2)) . 

Applying the Law of Sines to triangle AQB, we have 

sin(c) sin(p) 
sin(0) sin(/3) 

hence 

• f \ • (R\ ■ n • / ^ sin(p) sin(fl) . 

sm( a) sm(p J sm(c) = sm(a) — sm(c) 

sm(c) 

= sin (a) sin(p) sin(#) 

= 2 sin(a/2) cos(a/2) sin(p) sin(6>). 

The dual of Eriksson's result is [Tj 

t\ sin(a) sin(/3) sin(c) 



tan 
from which 



2/ cos(a) + cos(/3) + cos(7) — 1 



tan (L) = - ^W2)sin(p)sin((9) 

V2/ sin(a/2) + cos(p) sin(0) V ; 

follows immediately. 

Adding equation (JSJ) to (jlj, we obtain [9] 

cog( ) = com + cos( 7 ) _ cos (gg) cos (V) 
2sin(f)sin(0) sin (f ) sin(0) 

subtracting equation (jSJ) from (jlj), we obtain 

... cos(/3) - cos( 7 ) sin (gg) sin (^) 

cos(6>) = — = — / \ — • 

2cos(f) cos(f) 
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Thus, given a, (3, 7, it is easy to compute 9 and then p (in that order). 
Rearranging equation <^ to 

tan(r/2) sin(a/2) + (tan(r/2) cos(p) + cos(a/2) sin(p)) sin(0) = 0, 

that is, 

tan(r/2) sin(o;/2) 

sm(#) : 



tan(r/2) cos(p) + cos(a/2) sin(p) 

we obtain the expression for 9 = f(p). The smallest admissible value p = p t hrcs occurs 
when sin(#) = 1, that is, 

tan(r/2) cos(p) + cos(a/2) sin(p) + tan(r/2) sin(a/2) = 0. 

Solving 

y x + z\/l — x 2 + y\/l — z 2 = 

is made possible via 



.2 



y- (x+y/1- z 2 ) = z 2 (l-x 2 ), 
hence 

(y 2 + z 2 ) x 2 + 2y Vl - z 2 x + (y 2 - z 2 - y 2 z 2 ) = 

hence 

x = - — 



y 2 + z 2 1 + ^T^ 5 

llGIlCG 

. . cos(r/2) + sin(a/2) 
cos(p) — 



l + cos(r/2) sin(a/2) 
gives the desired threshold. 

7. Definite Integrals 
7.1. Crofton/Exhumatus. We wish to evaluate 



J 7r — arctan (J tan 2 | esc 2 | + 1 tan | J 

/ sin xdx if < y < ir, 

J ■ /tan 2 f esc 2 1 + 1 



7 arctan ( wtan 2 | esc 2 | + 1 tan 1 1 

/ —. = sin xdx if n < y < 2n 

1 . Ann 2 £ csc 2 I _u 1 



/tan 2 £ esc 2 Jf + 1 
A miraculous substitution 



COS Z = COS I COS |, < Z < 7T 
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is due to Crofton & Exhumatus [2]; from 



sin z dz = | sin | cos | dx 



we deduce that 



sin z dz cos z sin zdz 4 cos z sin z 

sm x da; = 2 sin 5 cos - y — — — — = 4 = — - dz 



2 sm | cos | cos | cos | cos z 2 



and 



tan 2 §csc 2 | + l taU " : 



tan 2 | 



because 

tan 2 z + 1 = sec 2 z = sec 2 | sec 2 | = (tan 2 | + lj sec 2 | 

= tan 2 § sec 2 | + (tan 2 \ + l) 

= (tan 2 f sec 2 | + tan 2 | ) + 1 

= (tan 2 § esc 2 | + tan 2 § + 1. 

Since cos(z), cos(y/2) obviously have the same sign, it follows that tan(z), tan(y/2) 
likewise have the same sign and 



tan z — \l tan 2 | esc 2 | + 1 tan |. 



If < y < 7r, clearly cos(?//2) > and the range < x < it maps to y/2 < z < 
ir/2. Also, tan(y/2) > 0, hence 



z = arctan ( ^/tan 2 | esc 2 | + 1 tan | 



hence 



7T - arctan (Jtan 2 § esc 2 § + 1 tan |) tang. ,4coszsinz 

■sinxax = -iix — z) — — dz 



tan 2 fcsc 2 f + l tanz cos 2 § 

4 tan | 



cos 2 I 



(7r — z) cos z dz. 



The required definite integral is thus 

7r/2 



4 tan | , 

— jj- I (it — z) cos zdz 

cos t ■/ 

2 y/2 
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which is elementary. 

If 7r < y < 27T, clearly cos(y/2) < and the range < x < it maps to y/2 > z > 
7r/2. Also, tan(y/2) < 0, hence 



z = 7r + arctan ( ^/tan 2 | esc 2 | + 1 tan | 



hence 



arctan ( wtan 2 | esc 2 | + 1 tan |J 

sin x <ic 

tan 2 | esc 2 | + 1 



tan | 
tan z 



(z-tt) 



4 cos z sin z 
cos 2 I 



4 tan | 2 

-(7T — cos z dz. 



COS 



2 ^ 



The required definite integral is thus identical to before (although here the lower limit 
y/2 is greater than the upper limit 7r/2). 
Computer algebra swiftly gives 



1 + 



d 
dy 



4 tan I f 2 

— j- I (7r-z)cos zdz 

cos^ - 



y/2 



(y 2 - Airy + 3ir 2 - 6) cos(y) - Q(y - 2vr) sin(y) - 2(y 2 - Arcy + 3vr 2 + 3) 

8cos(y/2) 4 

and this is useful at the conclusion of [3.1] . 

7.2. Jones/Benyon-Tinker. Combining our results with those in [U [5], we have 



r/2 



sin(r — k — p) sin(fi;) sin(p) 



t/Lk V sin ( K ) 2 sin(p) 2 - [cos(k) cos(p) - cos(r - k - p)\ 



-Ap (7) 



£;(sin(f))-cos(^) 2 x(sin(f 



2 / _ \2 

f()S(|J -C0S(^2" K 



sin(/t) 



in connection with primal perimeter |2.3] and 

7r-(re-o72) 



sin(cr — k, — 6 s ) sin(fi;) sin(#) 



cr/2 



sin(fi:) 2 sm(8) 2 — [cos(k) cos(#) — cos(cr — k 



--d6 (8) 



E ( cos (%)) -sinful if (cos 



sin (§) - sin 



sin k) 



Random Spherical Triangles 



27 



in connection with dual area |2.4j . A direct symbolic proof of these formulas is not 
known [3.2] . 

Consider the problem of integrating equation ([7j) with respect to k, < k < t/2 
and of integrating equation (|E|) with respect to k, a/2 < k < ir. In ((7j), p is integrated 
out first, k second. In ((H]), 9 is integrated out first, k second. By symmetry, we gain 
nothing by integrating out k first, thus a closed- form expression for unconditional 
density would seem unlikely. Another miraculous change of variables might, however, 
be brought into play. Other approaches based on other coordinate systems exist |5.1[ 
15.2] . It is still too early to rule out the possibility of a breakthrough here. 
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We are grateful to M. Larry Glasser for a helpful discussion about integrals at the 
end of |2.3j & |2.4j . Much more relevant material can be found at [TQ1 [EE], including 
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studies length distributions for open and closed random n-step tours {n > 3) on 
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